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The spectrum from Lattice NRQCD 

Presented by C. T. H. Davies a * 

a Dept. of Physics and Astronomy, University of Glasgow, Glasgow, G12 8QQ, Scotland 

I review recent results for heavy-heavy spectroscopy using Lattice NRQCD. The NRQCD collaboration reports 
that spin-independent splittings for the T are scaling for a sensible range of f3 values in the quenched approxima- 
tion. Spin-dependent splittings are not, if the scale is set by spin-independent splittings. Results which include 
higher order spin-dependent relativistic and discretisation corrections show differences from previous (NRQCD 
collaboration) results without these. As expected, the differences are small for T but rather large for charmonium. 
New results from the SESAM collaboration for T spectroscopy on configurations with Wilson dynamical fermions 
show good agreement with previous results on HEMCGC configurations with staggered dynamical fermions. 



1. New results 

The new results using NRQCD which I will dis- 
cuss are : 

• SESAM collaboration results for the T 
spectrum on configurations with 2 flavours 
of dynamical Wilson quarks. 

• Results from the NRQCD collaboration for 
the T spectrum on UKQCD quenched con- 
figurations at (3 — 6.2, giving a more de- 
tailed spectrum than previous results at this 
fine lattice spacing. 

• Results for the \P and T spectrum which 
include additional relativistic and discreti- 
sation corrections in spin-dependent terms 
from Trottier, SESAM and UKQCD. 

• Further B c results on dynamical configu- 
rations and including relativistic c quarks 
with non-relativistic 6s. 

2. NRQCD 

The splittings between radial and orbital ex- 
citations for systems made of heavy quarks are 
around 500 MeV, much less than the masses of 
the bound states. This implies that these are 
non-relativistic systems and a systematic expan- 
sion of the QCD Hamiltonian in powers of v 2 may 
be useful §. 

*Work supported by PPARC, NATO under grant 941259 
and EU network grant CHRX-CT92-0051 



The continuum action density, correct through 
0(Mqv 4 ), is broken down according to 
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HQ OUt and SH cont are given explicitly in ref.[Q. 
On the lattice Ho and the leading piece of 5H, 
SHi are given by: 
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The last two terms in bH\ come from finite lat- 
tice spacing corrections to the lattice Laplacian 
and lattice time derivative, of 0(a 2 MqV 4 ) and 
0(aMqV i ) respectively, n is the stability param- 
eter used in the evolution equation below. 

The quark propagators are determined from 
evolution equations that specify the propagator 
value, for t > 0, in terms of the value on the pre- 
vious timeslice; 
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The quark propagators are combined with 
smearing operators at source and sink to produce 
good overlap with different states. This is done 
in different ways by different groups, NRQCD 
and SESAM using Coulomb gauge wavefunction 
smearing and UKQCD, gauge-invariant blocking. 
All use multi-exponential fits to multiple cor- 
relation functions on large ensembles to obtain 
masses for radial excitations in s and p channels 
which can be compared to experiment. 
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Figure 1. The spin independent spectrum for bb 
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to be added to 8H\ above in the evolution equa- 
tion. The action used is then correct through 
0(Mqv 6 ) for spin- dependent terms. In addition 
these groups have included extra discretisation 
corrections for the spin-dependent terms in 5Hi 
at £>(M Q v 4 ), i.e. terms of 0(a 2 M^v 6 ). These 
involve |g{ replacing the E and B fields with an 
improved version: 
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and replacing the spatial derivative with an im- 
proved version (as was done for the leading spin- 
independent terms in 8H\) : 



A( 2 )=A( 2 )--A< 4 ). 
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These corrections then appear as additional terms 
in 5H2- The results discussed here differ in how 
groups have treated E and B fields and correc- 
tions to derivatives in the spin-independent terms 
at 0(Mqv a ). This should not cause a big ef- 
fect since these corrections are next-to-next-to- 
leading order in spin-independent splittings. 

The Ci and coefficients should be matched 
to full QCD either perturbatively or non- 
perturbatively. This has not been done by any 
of the groups. Instead they have relied upon 
tadpole-dominance arguments || to replace the 
gauge fields appearing in derivatives and E 
and B fields by U^/uq. uo represents the ef- 
fect of tadpoles in reducing the mean value of 
the link. Different values can be taken for uq 
and they will be compared below. After these 
modified gauge links are used in the action the Ci 
and fi are set their tree- level values of 1. There 
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Figure 2. T fine structure 



is evidence from pcrturbativc calculations of the 
NRQCD self-energy diagram that the radiative 
corrections to this, at least for c\ and C5, are 
very small (less than 10%) for reasonable values of 
Mqo- Q. However, they could in principle be of 
the same size as the relativistic and discretisation 
corrections of SH2- 

There are two parameters then to fix in the 
lattice calculation: a and Mq. For a, the 2S-1S 
or IP-IS splitting should be used, where S and 
P are spin averages over s and p states. This 
has the advantage of being independent of quark 
mass experimentally in the b, c region. Because 
not all the states have been seen experimentally 
the 2 3 S'i - l 3 _Si and 1 3 P - l 3 S"i splittings are 
used in the bb sector, where the 3 P is the spin 
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Table 1 

Parameters for the ensembles used in the results 
discussed. 



average of the 3 Po.i,2- To fix the quark masses the 
energy at finite momentum is calculated for one 
meson (say the l 3 ^) and the denominator of the 
kinetic term in the dispersion relation is taken as 
the absolute mass of that meson in lattice units. 
The energy at zero momentum differs from this 
since the mass term was dropped from Hq- The 
difference between the energy at zero momentum 
and the mass can be compared to perturbative 
predictions and agrees well || . 

3. Results for Upsilon spectroscopy 

The results to be discussed are tabulated below. 
Where a symbol is denoted in the last column 
that is the one used in Figures 1 and 2. 

The HEMCGC configurations use Kogut- 
Susskind dynamical fermions and have two en- 
sembles, one with ma = 0.01 and one with ma 
= 0.025. The NRQCD results on these configu- 
rations have been previously reported The 
matching quenched results at (3 — 6.0 represent 
a higher statistics study than that reported in 
ref. @. The SESAM configurations use Wil- 
son dynamical fermions and have three ensem- 
bles. The two with lightest dynamical mass (k = 
0.157 and 0.1575) were used for the T spectrum 
calculations by Achim Spitz and Henning Hoc- 
ber. SESAM results on these configurations are 
reported elsewhere in this Proceedings [0. 

Figure 1 shows the collected results for the 
spin-independent spectrum and Figure 2 the fine 
structure. The values of the lattice spacing used 
to convert to physical units have been taken as the 
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Figure 3. Thejatio of splittings 2 3 S'i - l 3 ^ to 
1 3 P — l 3 5i in bb as a function of lattice spacing in 
fm at 3 different values of the lattice spacing on 
quenched configurations (NRQCD collaboration). 
The solid line represents the experimental result. 



average of that from the 2S-1S and IP-IS split- 
tings. a -1 takes the value 2.4 GeV in all cases 
except for the heavier dynamical mass SESAM 
results in which it is 2.3 GeV. A bare b quark 
mass in lattice units of 1.71 was obtained by the 
NRQCD collaboration by tuning on the quenched 
[3 — 6.0 configurations. A value of 1.8 was used 
by them on the HEMCGC unquenched configu- 
rations but results from a kinetic mass analysis 
of the T H indicate that this was too large and 
1.7 would have been better. The SESAM results 
use a bare b mass of 1.7. 

Spin-independent spectrum - It is clear 
from the open circles in Figure 1 that the spin- 
independent spectrum on quenched configura- 
tions is not correct. We expect there to be er- 
rors because the coupling constant runs incor- 
rectly between the scales appropriate to, say, the 
IP and the 2S, so that it is not possible to fix an 
effective coupling which gives the correct answer 
for both states. Before comparing results at dif- 
ferent values of n/, however, we must first check 
that the results are scaling for a given value of 



Figure 4. The ratio of the 1 3 P - l 3 Si splitting 
in bb to the p mass as a function of lattice spacing 
in fm at 3 different values of the lattice spacing in 
the quenched approzimation. T results from the 
NRQCD collaboration. The circles use UKQCD 
results for the p mass, the squares GF11 results. 



n f . 

NRQCD is an effective theory reproducing the 
low energy behaviour of QCD but in which the 
ultra-violet cut-off plays a crucial role. It is there- 
fore not possible to take a to zero within NRQCD. 
The Ci and fi coefficients will start to diverge 
as powers of I/Mqci Q| and we will lose con- 
trol of the NRQCD expansion. However, there 
is no need to take a to zero if we can demon- 
strate a independence of our results for a reason- 
able range of values of a. Provided Mqcl > 0.6 
the Ci coefficients which have been calculated arc 
perfectly well behaved Q and we do not expect 
any problems. For the T system this corresponds 
to < 6.4. 

The NRQCD collaboration now has results at (3 
values of 5.7, 6.0 and 6.2 in the quenched approx- 
imation (see table 1) |J. Figure 3 shows the scal- 
ing of the ratio of the 2S-1S to IP-IS splittings. 
There is no sign of significant scaling violations 
and the result clearly disagrees with experiment. 
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Figure 5. The ratio of splittings 2 3 S'i - l 3 Si to 
1 3 P — l 3 £i in bb as a function of the number of 
dynamical flavours. NRQCD results are circles, 
SESAM results, a star. The solid line represents 
the experimental result. 



These results use the action with 8H\ described 
above in which the leading 0(Mqv 2 ) terms are 
corrected for their lowest order discretisation er- 
rors. The remaining leading discretisation error 
is then expected to be the difference of terms of 
0(a 4 M^v 6 ) between s and p states. At = 5.7 
the 0(a 2 ) errors arising from the gluon field con- 
figurations generated with the unimproved Wil- 
son plaquette action become significant. They 
can be corrected for perturbatively || and this 
has been done in Figure 3. It amounts to a 5% 
reduction in the IP- IS splittings at (3 — 5.7, less 
than one a in the ratio. 

Figure 4 shows the ratio of the IP- IS splitting 
in the T system to the UKQCD values of the p 
mass at the same three values of (3. The UKQCD 
results shown are those for the tadpole-improved 
clover [[l0| light fermions. Good scaling is seen 
and very clear disagreement with experiment in 
the quenched approximation. For comparison I 
also show the ratio using the GF11 unimproved 
light hadron results - clear violations of seal- 



Figure 6. The l 3 5i - l^o splitting in MeV with 
the scale set by the 2 3 S'i — l 3 Si splitting in bb as 
a function of the squared lattice spacing in fm 2 
at 3 different values of the lattice spacing in the 
quenched approximation. NRQCD collaboration; 
the open circles use uop as in the simulation, the 
open squares are rescaled to give the results using 
uql- The star gives UKQCD results with a higher 
order action and uqp. 



mg are seen. 

Having demonstrated scaling of the spin- 
independent spectrum at rif — (and therefore 
presumably at other values of n/ also), we can 
now study rif dependence of the results. Figure 
5 shows again the 2S-1S/1P-1S ratio plotted as a 
function of rif . The results at rif — 2 are in much 
closer agreement with experiment than at rif = 
0. The two results using different ensembles at 
rif = 2 with different types of dynamical quarks 
are in good agreement with each other. 

Since important momentum scales in T s and 
p states are around 1 GeV we might expect these 
splittings to 'see' 3 flavours of dynamical quarks. 
Extrapolating linearly through the NRQCD re- 
sults at rif = and 2 does cross the experimen- 
tal line at rif = 3. Improved statistics at rif = 
2 would be useful to show definitively that this 
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Figure 7. The l 3 5i - l^o splitting in MeV with 
the scale set by the 1 3 P2 — l 3 -Po splitting in bb as 
a function of the squared lattice spacing in fm 2 
at 3 different values of the lattice spacing in the 
quenched approximation. Open circles, NRQCD 
collaboration; star UKQCD collaboration with a 
higher order action. 



point gives incorrect answers. 

The numbers plotted at rif = 2 are for the light- 
est dynamical quark mass in both the NRQCD 
and SESAM cases. In principle the results should 
be extrapolated (possibly linearly ]1^]) in the 
light quark mass to the point 
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However, no significant md yn dependence has 
been seen in either dynamical ensemble and so 
this has not been attempted. 

Spin-dependent spectrum - Figure 2 com- 
pares the fine structure for the sets of results in 
Table 1, and includes details of the action used 
by each group. It is immediately obvious that the 
fine structure is much more sensitive to changes 
in the quark action. This is not surprising as 
the fine structure appears for the first time at 
an order one less in the relativistic expansion 
than the spin-independent spectrum. 5H\ used 



by the NRQCD collaboration contains only the 
leading order spin-dependent terms and so we ex- 
pect changes at the 10% level (v 2 « 0.1) on using 
5H2- It is also clear that unquenching has a big 
effect and will be necessary to get the right an- 
swers. 

Comparable results with and without 8H2 are 
the CO from UKQCD @ and the o from NRQCD, 
both on quenched configurations at p = 6.0. A 
significant (15er) effect is seen in the hyperfinc 
splitting and the shift corresponds to a 10% effect. 
The shifts in the p fine structure may turn out 
to be somewhat larger (possibly as much as 20%) 
but currently they are not statistically significant. 
A comparison of the two collaborations results 
with only 5H\ shows agreement between the two 
different methods of smearing. 

The txi points include both relativistic correc- 
tions (which are physical) and discretisation cor- 
rections (which are not). The scaling violations 
in the quenched NRQCD results can attempt to 
untangle these effects. Figure 6 shows a scaling 
plot of the hyperfine splitting divided by the 2S- 
1S splitting as a function of a 2 in fm 2 . We expect 
violations of scaling at this order in the NRQCD 
results from errors in the B field in the a ■ B term 
which gives rise to this splitting. Figure 6 shows 
clear scaling violations with a slope p, given by ~ 
900 MeV, when we write [111 : 
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25-15 
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It is not surprising to find scaling violations 
of this size since the hyperfine splitting responds 
to much shorter distances than spin-independent 
splittings and is given by a single term in SH 
without the cancellation that occurs for deriva- 
tive type discretisation terms in spin-independent 
splittings. 

The hyperfine splitting is also sensitive to the 
value of the bare quark mass (like X/Mqo) and 
errors in how well this is fixed will affect Figure 
6. 

In addition there is sensitivity to the value of 
Mo, because the B field contains 4 links. We ex- 
pect the hyperfine splitting to vary as 1/uq, and 
this was borne out by work in ref. || where results 
with uq and without (uq = 1) were compared. 
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NRQCD and UKQCD results both use uq from 
the fourth root of the plaquette, denoted uqp in 
Figure 2. Recent work by Trottier |H| has sug- 
gested that Mo from the Landau gauge link, uql, 
might provide a better estimate of the radiative 
corrections to C4. This might also be true for 
other a if uql captures tadpole effects more ac- 
curately. This value of uql is therefore used by 
the SESAM collaboration and has the effect of 
increasing the hyperfine splitting over the value 
that would be obtained with u 0P . 

Figure 6 gives NRQCD values for the ratio of 
hyperfine to 2S-1S splittings that might be ex- 
pected using uql by rescaling the results obtained 
with uqp by the sixth power of the ratio of uq 
values. The uql results are somewhat flatter and 
this might indicate that some of the previous scal- 
ing violations arose from radiative corrections to 
C4. Perturbative or non-perturbative calculations 
of various Cj will be needed to answer the question 
of which uq is better (if there is a single answer) 
and allow us to include radiative corrections in 
the leading coefficients in a consistent next-to- 
leading-order calculation of fine structure. 

Inn/ extrapolations of fine structure splittings 
it might not be true that rif = 3 is the relevant 
physical point, n/ = 4 is somewhat more likely 
given the short distance nature of these quanti- 
ties. In that case the extrapolations in n/ must 
be done in terms of other quantities for which n f 
= 4 is also the physical point, and not the 1P- 
1S or 2S-1S splittings. Indeed, there will be no 
value of degenerate dynamical rif for which we 
could get the right answer for the ratio of fine 
structure to spin-independent splittings. Instead, 
for example, we would have to extrapolate the 
ratio of the hyperfine splitting to one of the fine 
structure splittings. For this ratio better scaling 
is seen in the NRQCD results, as shown in Figure 
7, although the results at (3 — 6.2 are still rather 
uncertain. So an extrapolation in rif could then 
be done. 

The results at different rif are shown in Fig- 
ure 8, where NRQCD results are compared at 
rif = and 2, and UKQCD results at rif = 
are compared to SESAM results at rif = 2. The 
SESAM results are not strictly comparable with 
UKQCD since they use u l; if u p had been used 



Figure 8. The ratio of the l 3 ^ - l 1 ^ splitting 
to the l 3 i-2 — l 3 -fo splitting in bb as a function 
of rif. Open circles, NRQCD collaboration; star, 
UKQCD collaboration at rif =0, SESAM collab- 
oration at rif — 2. See text. 



their result for the ratio plotted would be smaller 
by about 5%. Unfortunately the fine structure 
splittings are not very accurate and the resulting 
error in extrapolating the hyperfine splitting in 
this way is very large. At rif =4, the hyperfine 
splitting would be 40(10) MeV. 

A comparison of the p fine structure to exper- 
iment from Figure 2 shows interesting features, 
albeit with large errors. The overall scale of the 
splittings given by 3 p2 — 3 Po is obviously too small 
in the quenched approximation (and is reduced 
further when scaling violations are removed as for 
the hyperfine splitting). The unquenched results 
are much closer to experiment. 

The ratio of splittings, ( 3 P 2 - 3 Pi)/( 3 Pi - 3 P ), 
is larger than experiment for the NRQCD results, 
1.1(4) versus 0.66(2). That this is probably a dis- 
cretisation error is borne out by the apparent im- 
provement in this ratio in the UKQCD results. 
The ratio viewed in a potential model picture has 
differences of long and short range effects in it and 
it would be unlikely to be correct in the quenched 
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Figure 9. The hyperfine splitting for charmo- 
nium versus a 2 taken from ref. |i"q |. The next- 
to-leading order NRQCD results are given with 
itop, open circles, and with uol, filled boxes. The 
leading order NRQCD result |l6| is marked with 
a cross. The relativistic heavy Wilson approach 
fl^l gives the open boxes. 



approximation, but is hard to determine accu- 
rately. The unquenched results in Figure 2 look 
encouragingly in agreement with experiment but 
a full study of scaling and n f extrapolations must 
be done. 

Obviously more work is required on the fine 
structure of the bb system and I believe that 
NRQCD will provide the most accurate results 
in this area. It will also be important to under- 
stand which quantities should be extrapolated to 
which values of n/ or, failing this, to make config- 
urations with 'real world' dynamical quark con- 
tent. Only then will we be able to get the correct 
answer for all the splittings without nf extrapo- 
lations. 

4. Results for Charmonium spectroscopy 

The charmonium sector is a difficult one to sim- 
ulate on the lattice, because it is neither very rel- 
ativistic or very non-relativistic. For NRQCD to 
remain above M c a of around 0.6 requires (3 < 5.85 
]16[| . Each order of relativistic correction is only 



30% smaller than the previous order since v 2 « 
0.3. For the standard relativistic heavy Wilson 
approach |i~7| the opposite requirement is true, 
(3 > 6.0, because of problems with fixing the me- 
son mass as a shift from the energy at zero mo- 
mentum |l8| l . 

There has been a long standing issue JllJ 
that results for charmonium hyperfine splittings 
disagreed between these two methods, with the 
NRQCD results giving a hyperfine splitting of 
around 100 MeV and the relativistic method 
around 70 MeV (experiment 116 MeV). 

New work by Trottier in NRQCD Q using the 
relativistic and discretisation corrections of 8H2 
has shown that indeed the relativistic corrections 
from higher order terms in NRQCD are large. He 
finds the hyperfine splitting becomes smaller be- 
tween the 8H1 results and the 8H2 results by a 
factor of 60%. This is certainly of O(30%) so 
not in principle surprising but nevertheless dis- 
appointing to proponents of NRQCD. 

A shift of the same size in spin-independent 
splittings would be a 15% effect (certainly of the 
same order as the naive expectation of 10%) since 
they are lower order in the non-relativistic expan- 
sion. This means a systematic error in a -1 deter- 
minations from NRQCD charmonium at this level 
if next-to-next-to-leading order spin-independent 
terms are not included. This is significantly less 
than the statistical error that can be achieved, 
but might nevertheless be acceptable in some ap- 
plications. 

It seems likely that future progress in the char- 
monium sector will use the heavy Wilson ap- 
proach, possibly corrected for p A terms, as has 
been suggested p0| . 

5. Results for B c spectroscopy 

Recent experimental evidence for the B c par- 
ticle encourages lattice predictions for the spec- 
trum of be bound states. Early work |2lJ has 
used quenched configurations and NRQCD for 
both the b and the c quarks. 

More recent results are collected in Figure 10. 
These include NRQCD collaboration results anal- 
ysed by Martin Gorbahn on dynamical configu- 
rations from the MILC collaboration ((3 = 5.415, 
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Figure 10. Lattice results for the spectrum of 
the B c system. Open circles are results using 
NRQCD for both b and c quarks on quenched 
configurations at (3 — 5.7. Filled circles use un- 
quenched MILC configurations. XI are results us- 
ing NRQCD b quarks and heavy Wilson c quarks 
on quenched configurations at /3 = 6.2. Error bars 
represent statistical uncertainties only. Dashed 
lines show results from a recent potential model 
calculation [E2fl 



rif = 2, KS, ma =0.0125) and UKQCD results 
from Hugh Shanahan using NRQCD b and rela- 
tivistic (tadpole-improved clover) c quarks on (3 
= 6.2 quenched configurations. For comparison 
the older results |^l| are shown and results from 
a potential model analysis |^2| . There is no clear 
discrepancy with the potential model results as 
yet, despite the fact that the c quark is more rel- 
ativistic in B c than in charmonium. 

One of the problems with a mixed system like 
the B c is an ambiguity in what to take for the 
quark masses. This is quite significant in the 
quenched approximation because the value of a -1 
depends on whether it is fixed from the T or ^ 
system and the values of the bare quark masses 
have also been fixed separately within these sys- 



tems. The splittings in Figure 10 will change 
somewhat if the quark masses were altered. To 
perform consistent extrapolations in n f it will be 
necessary to use one particular splitting in, say, 
the T system to fix a -1 and then fix m c using 
this a -1 in the system. This has not been done 
as yet. 

One interesting feature of the B c system is its 
similarity to heavy- light systems, allowing a test 
of some of the techniques that will be useful for 
the spectrum there. In particular the spin 1 p 
states will mix because of a lack of charge conju- 
gation. It was possible for the quenched NRQCD 
results [^lj to diagonalise the mixing matrix and 
pick out the physical 1 + and 1 + states. This will 
also need to be done in the B sector, but has not 
been possible as yet p3[ . 
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